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1. INTRODUCTION 

Energy is a major factor in the growth of industrial production. Renewable energies mainly wind 
and solar are becoming the most supported solutions; as a result of less abundant fossil fuels, greenhouse gas 
effect and high population growth [1]. On the other side, wind power is a globally available solution unlike 
solar power due to its limitation to specific regions. Wind energy conversion systems are the fastest expanding 
technology and the key components to clean and save energy [2]. Among all variable-speed wind turbine 
designs, two commercial leading technologies are introduced in extensive scale; wind turbine based DFIG and 
synchronous generator (SG). Indeed, DFIG technology has overtaken other ones because of its cost, robustness, 
reliability and small power converters requirements [3], [4]. 

In this work, we are concerned with the control of the whole wind energy conversion system formed 
by a wind turbine connected to the grid via a DFIG, a back-to-back converter and three phase RL filter. The 
back-to-back converter includes two identical two-level voltage source converters; a rotor side converter (RSC) 
and a grid side converter (GSC), linked by a capacitor. The converter consists of six insulated gate bipolar 
transistors (IGBT) with anti-parallel diodes to allow the flow of current in both directions [5], [6]. To get full 
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benefits from the DFIG machine, different control techniques are employed for the synthesis of the DFIG wind 
turbine system such as the field-oriented control used to decouple and regulate the torque and the stator reactive 
power with classical PI-regulator, proposed in [7], [8]. DTC/DPC strategies are employed to control the WECS 
powers in [9], a fuzzy logic controller is designed for power capture optimization in [10]. To optimize power 
production of DFIG, a host of research works controlled the discussed system using nonlinear control strategies 
such as the backstepping approach designed in [11], feedback linearization technique presented in [12], and 
flatness-based approach used in [13]. However, in these works, the mechanical torque and the generator speed 
are considered as measurable variables which is inappropriate in industrial solutions for electric drive systems 
as long as one seeks to eliminate the mechanical sensors because they increase the maintenance costs and 
reduce the system reliability, thus affecting the overall wind turbine system stability. In [14], authors proposed 
a sensorless backstepping control technique for the DFIG. Nevertheless, only the generator speed is estimated 
whereas the mechanical torque is assumed to be measurable. In this paper a high gain observer is synthesized, 
in order to provide online estimation of the mechanical torque and the generator speed, only the electrical 
variables are assumed to be as measurable quantities. 

This article presents a sliding mode controller design in order to control the power converters associ- 
ated with the generator which are the key components to handle the energy between the wind turbine and the 
electrical grid. The RSC aims at regulating the generator speed to track the optimal reference provided by the 
MPPT block, the speed controller is designed based on the mechanical observer, and regulate the stator reactive 
power so as to be null. The GSC aims at guaranteeing that the DC-link is maintained at a given constant value 
and ensuring a unitary power factor by keeping the grid side reactive power at zero. 

The remaining sections of the paper are organized as follows: In Section 2, the whole wind energy 
conversion system mathematical model is presented, followed by the sliding mode controller design described 
in details for the proposed system in Section 3. The high gain observer is introduced in Section 4. The 
simulation results are discussed in Section 5. Lastly, conclusions and future scopes are drawn in Section 6. 
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Figure 1. Wind energy conversion system topology 


2. SYSTEM MODELLING 

The complete system scheme is depicted in Figure 1. According to [5], [15], the mathematical model 
in state space form in d-q reference frame of the DFIG connected to the back-to-back converter can be summa- 
rized by the following differential equations: 


Tı = eb ee S —C, (1) 
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te = Ta (Žo = Rote q= wsLgzs) (6) 
g 
where |£1, £2, £3, £4, £5, £6] = [Om, irq, ird, Vdc, İgq, İga] denotes the state vector. wg, Uug, Va and vg are the 


back-to-back converter control inputs in d-q reference frame. Qm, Cg, p, J and fy are the DFIG mechanical 
rotational speed, mechanical torque, the number of pole pairs, the rotor inertia and the friction coefficient, 
respectively. (with Q+ = Qm/N and C; = C.N, N is the gear box gain. Since the shaft rotor coupling 
involves a gearbox). Ws, Wm and w, denote the stator flux synchronous speed, the angular rotor speed and 
the rotor angular frequency, respectively, (with wy = ws — Wm). ira, tga are the rotor and filter direct current 
components. trg, tgq are the rotor and filter quadrature current components. V, is the magnetizing amplitude. 
Rs, Rr and Rg designate the stator, rotor and filter resistor, respectively. L,,L,, Ll and Lg indicate the 
stator, rotor, mutual and filter inductance, respectively. Cg, is the DC-link capacitor, irae represents the current 
flowing through the RSC. u = (—3pLmVs)/(2Lsws) designates the electromagnetic torque constant and 
o = 1 — LŽ /(LsLr) is Blondel’s dispersion coefficient. 


3. CONTROLLER DESIGN 
3.1. Rotor side converter controller 

The RSC employs the sliding mode control technique to transfer the available electrical power from 
the DFIG stator to the utility grid. The electrical power is generated via the mechanical power produced by the 
wind turbine rotor. The controller involves two control loops: the first one is the speed control loop that aims 
at regulating the generator speed to track the optimal reference provided by the MPPT block. The second one 
is the stator reactive power control loop which is intended to control the stator reactive power so as to be null 
[16]-[18]. 


3.1.1. Generator speed control loop: 
Define the following second order sliding mode surface [19], [20]: 
Salt Cg) = Ceir óe (7) 


where e} = xı — xj denotes the tracking error between the generator speed and its optimal value, and 0, is 
the positive parameter design. 
It follows from (1) that the error e; undergoes the following equation: 
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Based on (1), (2), (7) and (8), the time derivative of the sliding surface is given by: 
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By setting the dynamics of the sliding surface to zero, one can obtain the equivalent control law which ensures 
the attractiveness and invariance of the sliding surface S,,: 
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Let us consider the positive Lyapunov function candidate: 


Wı = x (12) 
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Differentiating (12), yields: 
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By substituting (11) in (13), one gets: 
. UX4 
= w n 14 


The following discrete control law will be chosen to guarantee the dynamics of the candidate Lyapunov func- 


tion to be negative: 
20d Ly 
T af kusgn (Su) (15) 
[XA 
where ką is the controller parameter to be designed. 
Then, the generator speed global control law is given by: 
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3.1.2. Stator reactive power control loop: 
To design a controller for the stator reactive power exchanged with the grid, the sliding surface is 
defined as follows: 











Solz) = e2 = Q: — Q% (17) 
where the reactive power equation is given by [5]: 
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Using equations (3), (17) and (18), the sliding surface Sg undergoes the differential equation: 
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with So = 0, the equivalent control law Ud.eg Can be expressed by: 
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The following discrete control law will be constructed which ensures the negativity of the derivative of the 
candidate Lyapunov function W3: 

4 oL Ll, 
o 3 LimVsXa4 
where kg is the controller parameter to be designed. 
Consequently, the final control law for the stator reactive power is given by: 


4 oL.L 
= — (Rr — WO Lp = 
Ud s (Rr£3 — wro Lre) + ARRS 


Uim kosgn (SQ) (24) 


(kosgn (SQ) — Q:) (25) 


Int J Pow Elec & Dri Syst, Vol. 12, No. 2, June 2021 : 1139 — 1149 


Int J Pow Elec & Dri Syst ISSN: 2088-8694 g 1143 


3.2. Grid side converter controller 


The GSC also employs the sliding mode control technique, in order to exchange with the grid only 
the active power extracted or injected by the RSC from the rotor, which is indirectly controlled through the 
DC-link regulator [21], [22]. To this end, the GSC controller involves two control loops. DC-link and d-axis 
current control loops. On one hand, the DC-link control loop aims at regulating the DC-link voltage at a given 
constant value which is equal to the converter nominal voltage to avoid the damage of the RSC. On the other 
hand, the d-axis control loop is used to force the reactive power to be null, thus ensuring a unitary power factor 
by maintaining the direct current at its reference value that is equal to zero. 


3.2.1. DC-link voltage control loop: 
Let us define the DC-link voltage tracking error between the actual voltage and its reference as: 


e3 = £4 — Tå (26) 


Using (4), the dynamics of the error e3 are written as: 


a ee ee 3 Vs l r 
OSS (~5 a: F irae) — t4 (27) 


Taking into account the tracking error (26) and its dynamics (27), the considered sliding surface is defined: 
Ss (x) = €3 + û2€3 (28) 


where 02 is the positive parameter design. 
It follows from (5), (27) and (28) that the time derivative of the sliding surface S, has the next form: 
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In order to find out the equivalent control law that ensures the attractiveness and invariance of the sliding 
surface, we set the dynamics of the sliding surface S, to zero. 


AlO 
Ug.eq = a ap (x,t) (31) 


The Lyapunov function candidate can be defined as W3 = 1S? and in order to make its derivative negative, 
the following discrete control law will be constructed: 


ERON 
Ta 
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where k, is the controller parameter. 
In view of (31) and (32), the global control law for the DC-link voltage is given by: 
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3.2.2. D-axis current control loop: 

To control the grid side d-axis current 7g, we proceed in the same way as for the earlier control 
loops. Based on the alignment discussed in [5], one can calculate the reference of the d-axis current such as 
xg = (2Q7)/(3V;), where Q% is the grid side reactive power reference. 

Let us consider the sliding surface S;q(~) = e4 = x6 — x§. Using (6) the derivative of the sliding surface with 
respect to time is given by: 


, 1 
C=. (ue — Ryze + ws Ligits — a (34) 


Setting the sliding surface Sia = 0, one obtains the equivalent control law as: 


2 Ibas 
Vang a (Rite — Wslg%t5) + =, 8 (35) 


In order to get the discrete control law, the chosen Lyapunov function is defined as W; = 552). Using (34), 
the dynamics of the Lyapunov candidate are expressed as: 
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It can easily be checked that the discrete control law vgn = eee kiasgn (Sia) guarantees the negativity of 


Wa, where kia is the controller parameter. 
Therefore, the global control law for d-axis control loop is written as: 


2 2A T pig 
v == (Rete U: Lgts) F — (x6 — kiasgn (Sia)) (37) 
£4 T4 
Remark: The chattering phenomenon can occur by using the signum function sgn(.). In this paper, we suggest 
modifying the control laws by means of using a sigmoid tanh(.) function instead of sgn(.) in order to avoid 
bringing in the vibration and nonlinearity into control signals during the reaching mode. 


4. HIGH GAIN OBSERVER 

In this section, a mechanical observer is employed to get online estimation of the mechanical torque 
and the generator speed so as to decrease the number of sensors that reduce the reliability of the system. Since 
the DFIG’s dynamic behavior belongs to a class of relatively fast systems, the high gain observer is among the 
most suitable choices in the problem of state estimation. This choice is motivated by the fact that the observer 
structure involves an explicit correction gain [23]. To this end, a high gain state observer will be designed. Due 
to the fact that the position of the rotor is considered to be unknown because it needs a mechanical sensor, the 
d-q frame cannot take place as it needs the rotor position information. Thus, the compact form of the DFIG 
model in the a-6 reference frame, is constructed as in (38) and (39). 


w=Aw+g(u,i,u) +4 (t) (38) 
where w = (wy we w3) = Ce Cin. VE G) r is the state space vector, y is the output vector, 
0 fillies) 0 f2 (tas, Uap) P 
A=|0 0 -4),9(w,i,u) =| 4u,-2uw.],6@)=(0 0 e(t)),C=(1 0 0). 
0 0 0 0 


; ; ; ; : ; ; ; T T 
with w = (3/4) pLmis,gTotas: tag = (isa ls8 lra inp) » Uab = es UsB Ura urg) > 


fi (tap) = (30 ly) / 26 Lebi Lo tap: fo (Tap; Uag) = (SPLm Caleb, To (Fouag a Fiag), 
8 8 


p (whiz —LrLbmde\ p (Erh Lbmk\ p (—Rsbrl2  Rrbmle 
EN e a eded = tcl 2 Rite ds. RG) 


Int J Pow Elec & Dri Syst, Vol. 12, No. 2, June 2021 : 1139 — 1149 


Int J Pow Elec & Dri Syst ISSN: 2088-8694 g 1145 


[Q R [1 0  /Q -~i  /Q0 0 
n=(3, aezk t= ( o)l o) 
And £ (t) is a bounded function. Then, the state observer of systems (33) is designed as follows [24]: 


Ü = A û +g(Û, i u) — T (i) A, KC (ù -— w) (40) 
where T (i) = diag (1 fi (iag) —4fı (iag)) and Ao = diag (1 4 jz), 4 isa constant positive gain. 
The gain matrix K is selected so that the matrix (A — KC) is Hurwitz, with A = T7 +AT, which implies that 
P(A —KC)+(A -KOY P=-h, PER’. 

The high gain observer stability can be proved using the Lyapunov stability theorem and by considering the 
quadratic Lyapunov function candidate: Wy = £t P e, where £ =T (i) Apw and w = ù — w [24]. 


5. SIMULATION RESULTS 

The focus in this work is made on the design of an output feedback controller that is carried out 
using the nonlinear model of the entire wind energy conversion chain adopted. The proposed controller is 
developed based on the sliding mode technique and using a high gain observer. This latter is employed to 
estimate the generator speed and the mechanical torque, which improve the WECS performance and reliability 
via eliminating the mechanical sensors. In this study, we are interested in the control of the back-to-back 
converter, 1.e., RSC and GSC, taking into account the IGBT switching devices which are considered ideal 
in several works. In addition, our main concern in this paper is directed to both RSC and GSC making our 
study complete in comparison to other works that focus only on one of the two converters. Actually, the 
power transfer through the back-to-back converter can affect and disturb the global power injected into the 
utility grid. The model of the 2-MW DFIG-based WECS connected to the electrical network, employing the 
proposed output-feedback controller which was discussed in the previous sections and described by Figure 2 is 
performed using MATLAB/SIMULINK/SimPowerSystems environment. 
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Figure 2. Schematic representation of the proposed wind energy conversion system with the controller design 


The system parameters are given as follows: 

Wind turbine parameters: N = 100, R = 42m, p = 1.1225K g.m’. DFIG and filter parameters: Vp /Uņn = 
400/690V, f = 50Hz, Pa = 2MW, R,/R, = 2.6mQ/2.9mQ, Ls/Lr = 2.6mH/2.6mH, Lm = 2.5mH, 
R, = 0.02mQ, Ly = 0.4mH, C = 80mF. Observer and controller parameters: K = |25, 8, 33N, 0 = 90, 
ku = 250, kg = 300, ky = 120, kia = 150, 01 = 100, 02 = 100. 

To demonstrate the validity of the output feedback controller, simulation results are obtained in the presence of 
a realistic wind speed profile with an average wind speed of 9 m/s characterized by a 4% turbulence intensity 
over the time interval [0 100 s], as shown in Figure 3 (a), to further prove the efficiency of the proposed 
controller. Figures 3 (b)-(c), illustrate the generator speed and the mechanical torque of the DFIG and highlight 
the output feedback control loop convergence. Particularly, Figure 3 (b) shows the time-variations of the 
generator speed reference delivered by the MPPT block, together with the measured generator speed and its 
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estimated value. As it can be noticed, the estimated value first matches the measured generator speed and 
both track the reference effectively so as to meet the control objectives. Also, Figure 3 (c) displays that the 
mechanical torque value provided by the high gain observer meets the measured mechanical torque which 
demonstrates the efficiency of the mechanical observer. Figure 3 (d) proves that the stator reactive power 
converges to the suggested reference which implies that only the active power is exchanged with the grid on 
the stator-grid connection. This keeps the currents sinusoidal and in opposite phase with the grid voltages, as 
shown in Figure 3 (g), complying with the PFC requirements. Figure 3 (e) verifies that the DC-link voltage is 
maintained at its recommended level, which guarantees the appropriate operation of the converter on the grid 
side. Figure 3 (f) shows the convergence of the grid side d-axis current which ensures that the reactive power 
exchanged with the electrical network on the rotor-grid connection is zero, thus proving the achievement of 
the GSC objectives. Figure 3 (g) represents one phase of the network voltage and current. The latter depends 
on the variation of the generated active power which fluctuates according to the wind velocity. Figure 3 (h) 
illustrates the mechanical power extracted by the wind turbine blades and the electrical power exchanged with 
the utility grid. 
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Figure 3. Wind energy conversion system estimation and convergence 


The observer is synthesized using the DFIG model in the stationary reference frame (a-( frame) 
where the rotor position is not used for the observer modelling and implementation in order to increase the 
reliability of the system. The estimated value of the generator speed and mechanical torque are used as inputs 
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of the sliding mode controller that aims at improving the performance of the conversion chain so as to meet the 
control objectives. The used controller ensures the system to converge faster than other regulators [7], [8], [10]. 
In addition, it allows the system to work without the presence of sensors, minimizing the cost and complexity 
of the overall structure, unlike other works that must lean on the signals provided by the sensors to achieve 
their control objectives [15], [25]. Another added-value of the present study relies on extracting the maximum 
available power given by the wind turbine blades. It is worth mentioning that the reactive power ripples shown 
in Figure 3 (d) according to the proposed design are much smaller than those shown in [10], [16]. In fact, the 
ripples in active and reactive powers have an influence on the current injected into the electrical grid and on the 
stress on the turbine shaft. This can lead to the destruction of the wind turbine. 


6. CONCLUSION 


In this article, a sliding mode controller is designed for the complete WECS constituted by a wind 
turbine connected to the grid via a DFIG and a back-to-back converter. The advantage of the proposed design 
is that it eliminates the requirement of mechanical sensors using a high gain observer, which improves the 
efficiency of the controlled system. It has been illustrated by a numerical simulation that the performance of 
the WECS based on the sensorless-sliding mode controller, under realistic wind profile, is better compared to 
other controllers and it can be concluded that the control objectives are efficiently achieved. We look forward 
to continuing this study in our next papers, first, by developing a novel MPPT strategy to improve the power 
extraction without the need of using an anemometer sensor for wind speed measurement, and then, by testing 
the control laws in real-time simulations in order to validate the system’s performance. 
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